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We use a power-series expansion to calculate the eigenvalues of anharmonic os- 
cillators bounded by two infinite walls. We show that for large finite values of the 
separation of the walls, the calculated eigenvalues are of the same high accuracy as 
the values recently obtained for the unbounded case by the inner-product quanti- 
zation method. We also apply our method to the Morse potential. The eigenvalues 
obtained in this case are in excellent agreement with the exact values for the un- 
bounded Morse potential. 
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1 Introduction 

The spectra of many important potential functions frequently encountered in quantum 
mechanics can not be obtained exactly. Moreover, in most cases, the conventional ap- 
proximate methods commonly discussed in most standard textbooks(see, for example, 
[1]), while possibly leading to rough estimates, are either unsatisfactory or computation- 
ally complicated. A well-known important example in this respect is a sufficiently deep 
double-well potential, in which splitting of energy levels of two neighboring states, in par- 
ticular the ground and first excited state, is extremely small because of the non-degenerate 
nature of energy eigenvalues of bound states for one-dimensional potentials [1, 2]. Per- 
turbation and semiclassical methods are of limited use in such a case, and the variational 
method, although providing upper bounds, requires carefully chosen trial functions. A 
need thus arises for a relatively simple and effective approximate method with a high 
degree of accuracy. 

A variant of approximate methods and numerical techniques has recently been devised 
to calculate to high precision the spectrum of the one-dimensional symmetric anharmonic 
oscillators, with either a symmetric solution [3]-[7] or a symmetry-broken solution [8]-[10]. 



A method based on the Hill determinant has been developed in a recent study [8]. 
This method, known as the inner-product quantization procedure, uses the representation 
^(x) = J2i a i(E) x l Rp{x) of the wave function, where Rp is an appropriate reference 
function. With a heuristic argument based on the Hill determinant method, it is shown 
that the zeros of the coefficient functions ai{E) approximate the exact bound state energies 
with increasing accuracy as % — > oo. The method has been applied to several unbounded 
one-dimensional anharmonic oscillators, leading to highly accurate energy eigenvalues. In 
particular, for a sufficiently deep double-well, the splitting of the energies of the ground 
state energy and first excited state becomes apparent only after 26 significant digits. 

It has already, however, been pointed out [4, 5] that not all eigenvalues obtained 
by the Hill determinant method should be allowed, since the boundary condition that 
— > as \x\ — > oo is not incorporated into the method and may lead to incorrect 
results for some values of the coupling constants and for potentials in which x = oo is 
an irregular singular point of the Schrodinger equation. In addition, the choice of the 
arbitrary reference function is important, as emphasized in [8]: in general, it should fall 
off more slowly than the asymptotic form of the wave function. 

The purpose of the present work is to apply the method of power-series expansion 
to some anharmonic potential functions bounded by two infinite walls. We show that 
the expansion of the wave function in the form ty(x) = J2i a i x% in the finite interval 
— L < x < L with infinite walls at x — ±L leads, even for moderate values of L, to 
the same high accuracy as has been achieved for the anharmonic potentials considered 
in [8]. The initial work in this direction was carried out by Barakat and Rosner [3], 
for the case of a pure quartic oscillator. These investigators showed that the lower- 
order eigenvalues tend rapidly to the values of the unbounded oscillator as L is made 
large. The method has subsequently been employed for the bounded pure anharmonic 
oscillators x 2k , with k = l,---,5, and for shallow doubly anharmonic oscillators [4]. 
The power series in the finite range does not require introducing reference functions and 
accommodates all eigenvalues, since the boundary conditions are imposed at finite x. The 
problem concerning the convergence of the wave functions mentioned in [8] does not arise 
here. 

The rest of the paper is organized as follows. In Sect. 2, we apply the method of 
power-series expansion to anharmonic potentials, specifically to a quartic, sextic, octic, 
dectic, duodectic and quartic double- well. Also in Sect. 2, we consider potential functions 
admitting power series-expansion, as exemplified by the Morse potential. In Sect. 3, 
we explain the method for obtaining the energy eigenvalues numerically, presenting our 
eigenvalue calculations and comparing them with the values obtained by the inner-product 
quantization method and with the exact values available for the Morse potential. In 
Sect. 4, we present our conclusions. 



2 Power-series expansion solution 



To calculate the eigenvalues and eigenfunctions of the anharmonic oscillators bounded by 
infinitely high potentials at x — ±L, it is necessary to solve the eigenvalue equation (we 
here assume units K = 1 , 2 m = 1) 

^ + E-V(x) *(x) = 0, (1) 



with the boundary conditions \I/(±L) = 0. In the present work, we consider polynomial 
potential functions and potentials admitting power-series expansion. In the former case, 
our potentials V(x) are of the form 



V(x)=^x 2 + gx 2k , (A; = 2,.-., 6) (2) 

Here the coupling constant g > and the mass parameter /i 2 takes real values. 
Using the power-series expansion 



OO 



* = £a^ (3) 

1=0 

in eq. (1), we obtain the following recurrence formula for the expansion coefficients: 

ai = TIT^T) ' ^ M 

ai = 0, I < 

(4) 

The symmetry of eq. (1) implies that the solutions fall into two classes, even and odd. 
Even solutions can be obtained by imposing ao = l,ai = 0, and odd solutions can 
be obtained by imposing a = 0, a\ = l.(We ignore normalization in both cases.) 
The energy eigenvalues E are obtained from the condition ^(L) = for both even and 
odd solutions. Since we are dealing with potentials admitting power-series expansion for 
\x\ < L, the power-series solutions of ty(x) are convergent, by a well-known theorem in 
differential equations [11]. 

To illustrate the method for potential functions admitting power-series expansion, we 
choose the Morse potential. For this potential, which happens to be non symmetric, we 
can derive recurrence relations in a manner similar to the symmetric case(performing 
a power-series expansion for the potential and truncating it after certain power). The 
resulting power-series solution depends linearly on both ao and a\ and in addition depends 
on E. Extracting the coefficient of ao an d a,\ in the power-series solution and calling them 
f° and f 1 , respectively, we then cast the wave function *&(x, E) into a form suitable for 
applying the boundary conditions: 

E) = f°(x, E) a + f\x, E) a x . (5) 

The boundary conditions for \P(±L, E) yield two homogeneous linear equations with a 
and ai as unknowns: 

f°(L,E)a + f 1 (L,E)a 1 = 
f°(-L,E)a + f 1 (-L,E)a 1 = 0. 

(6) 

For there to exist a nontrivial solution for {a ,ai}, the following condition must be sat- 
isfied: 



f(L, E) f\-L, E) - f\L, E) f°(-L, E) = 0. 



(7) 



This fact allows us to solve for the unknown energies E and the constants a$ and a\. For 
the case of the well-known Morse potential 



V(x) = V (l-e- Xx ) 2 , 



(8) 



where V and A are the depth and range parameters, respectively, the exact energy eigen- 
values and eigenf unctions are given by 



E n = 2XJV 
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where n must be positive and less than — |, and where 



(9) 



y = 2ae 



-Xx 



The functions L a n (x) are the generalized Laguerre polynomials, while the normalization 
constant N(n) is explicitly given by 

"(»> = y^^iiF- 

To employ the method, we expand the potential in eq. 8 as power series, truncating 
the series beyond, say, x 30 : 
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3 Calculation of Energy Eigenvalues 

To obtain the energy eigenvalues, we approximate the power series with a finite number 
of terms ^ n (x, E). The boundary condition for a specific value of L is ^ n (L, E) = 0. To 
find the zeros of fy n (L,E) with respect to E, we first plot ^f n (L,E) as a function of E. 
Around each zero, we locate two nearby points E at which the function ^ n (L, E) changes 
sign. We use these two points in the initial iteration of the "bisection method" for finding 
the zeros. In implementing the method, we use the Mathematica (version 3) package, 
relying extensively on its exact-number manipulation capabilities. We check the stability 
of the numerical results at a certain degree of accuracy for a particular L by increasing n 
until the obtained value of E stays fixed. 

To illustrate the method, we present in Fig. 1 an example for the behaviour of even 
and odd solutions ip e {E) and ip {E) for L = 2 and I = 100 (where I indicates the number 
of nonvanishing terms in the truncated expansion of the wave function), for the case of 
the quartic anharmonic oscillator (potential V(x) = x 2 + x 4 ). 

We use the same procedure to obtain energies in the case of a non symmetric potential. 
In this case, however, we apply the condition given by eq. 7. 



In Table 1, we present the calculated energies of the ground and first excited states for 
the bounded quartic anharmonic oscillator. Our method shows systematic convergence for 
increasing values of both I and L, exceeding some of the high- accuracy results in [8]. As 
providing evidence for the accuracy attained, we present the ground-state energy with 300 
digits in Table 1. We have checked the rate of convergence with a reference function of the 
type used in ref. [8]. Unfortunately, however, our study reveals that this reference-function 
procedure does not accelerate the rate of convergence, and that it makes the determination 
of the energy eigenvalues less reliable. Table 2 shows the calculated energies of the ground 
and first excited states for the bounded sextic, octic, dectic and doudectic anharmonic 
potentials, in a comparison with the values reported in [8] for the unbounded case. 



Figure 1: The even solution %p e (E) (broken curve) and odd solution ip Q (E) (solid curve) for the bounded 
quartic anharomnic oscillator (potential V(x) = x 2 + x 4 ), where 2L is the separation of the walls and / 
is the number of non vanishing terms in the truncated expansion of the wave function. 
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Figure 2 shows the dependence of the ground state energy for the bounded quartic 
anharmonic oscillator (V(x) = x 2 + g x 4 ) on the coupling constant g, for < g < 10, 
and L = 5. It is clear that the ground-state energy increases, as expected, when the 
coupling constant is increased. 

Figure 3 shows the wave functions for the ground state and first excited state in the 
case of the quartic anharmonic oscillator (V(x) = x 2 + gx 4 ) for g — 0, 1, 2. Figure 4 shows 
the ground and first excited state for the sextic, octic, dectic and duodectic anharmonic 
potentials V(x) = x 2 + x 2k (where k — 3, 4, 5, 6, respectively). 

As emphasized in Sect. 2, all these graphs of the wave functions show the convergence 
of the power-series expansion in the finite interval. 

Table 3 presents the calculated energies of the ground and first excited states for 
the bounded double-well potential, in comparison with the values reported in [8] for the 
unbounded case. Our calculations confirm the predictions made in [8] concerning the 
pseudo-degenerate nature of the ground-state and first-excited-state energies for fi 2 = 25, 
which becomes apparent only after 26 significant digits. In addition, we present the 
energies of the ground state and first excited state for /i 2 = 35, 50. The splitting for 
(j 2 = 35 and fj 2 = 50 are revealed after 42 and 72 significant digits, respectively. 

Our calculations confirm the general qualitative results that the energies of the low- 
lying states of the double-well potential become almost degenerate (pseudo-degenerate) 
as the depth of the well increases. 

In the case of the Morse potential with the particular values A = 1 and V = 400, 



Table 1: The calculated energies for the ground and first excited state for the bounded quartic an- 
harmonic (potential V(x) = x 2 + x 4 ), where (2L) is the width of the well and I is the number of the 
nonvanishing terms in the truncated series of the wave function. 
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Table 2: The calculated energies for the ground and first excited state for the bounded quartic, sextic, 
octic, dectic and duodectic anharmonic oscillators (with the indicated potentials V(x)), where 2L is the 
width of the well, and where the underlined values are as calculated for the unbounded potential in [8]. 
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Figure 2: Dependence of the ground state energy of the bounded quartic anharmonic oscillator (potential 
V(x) = x 2 + gx 4 ) on the coupling constant. 




Figure 3: The wave functions for the ground (left) and first excited (right) for selected values of the 
coupling constant g, for the bounded quartic anharmonic oscillator (potential V(x) = x 2 + gx A ). 
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Figure 4: The solid curve represent the ground (left) and first excited (right) state wave function for 
the bounded sextic anharmonic oscillator (V(x) = x 2 + x 6 ). The dotted curves represent respectively 
(inward) the ground state wave functions for the bounded octic, dectic and duodectic case (V(x) = 
x 2 +x 2k , k = 4,5,6). 
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Table 3: The calculated energies for the ground and first excited state for the bounded potential 
V{x) — ~[i 2 x 2 + x 4 , where fi is the mass parameter and 2L is the width of the well, and where the 
underlined values are as calculated for the unbounded potential in [8] . 
/i 2 Parity E± 

L = 6 

+ 1.06036209 04841828 99647046 016 69266 35455152 08728529 

3.79967302 98013941 68783094 188 51256 89577660 65467327 

1 + 0.65765300 51807151 23059021 723 11059 35603749 37941936 

2.83453620 21193042 14654676 208 74896 49582169 40153453 
5 + -3.41014276 12398294 75297709 653 52190 91987123 39047565 
-3.25067536 22892359 80228513 775 54773 68771546 01147639 

L = 8 

10 + -20.63357670 29477991 49958554 83743150 87653159 46057736 
-20.63354688 44049110 79343874 10046139 03678429 34101495 



15 + -50.84138728 43819543 66250996 515 74123 37747896 27482985 
-50.84138728 41870051 54710149 735 64863 44459057 68683578 

25 + -149.21945614 21908880 29163966 538 16577 44754406 92275913 
-149.21945614 21908880 29163958 97435901 91957349 04923409 



35 + -297.91219449 37076184 50107115 88417483 94538893 52429871 
-297.91219449 37076184 50107115 88417483 94538891 68865338 

L = 9 

50 + -615.02009090 27578165 66217383 21036156 72635810 55195838 
06693413 44260555 97217697 29429892 
-615.02009090 27578165 66217383 21036156 72635810 55195838 
06693413 44260555 97217403 20140809 



we obtain the energy eigenvalues tabulated in Table 4. Our values are identical with 
the exact values to the order indicated in the table. Figure. 5 shows the ground state 
and the first excited state, respectively, for a bounded Morse potential confined to the 
interval — 2 < x < 2. These figures match the exact corresponding wave function of the 
unbounded Morse potential in the same interval. 

Table 4: Energy eigenvalues for the Morse potential with depth parameter Vb = 400 and range parameter 
A = 1, where the results for the bounded case are obtained for / (the number of nonvanishing terms in 
the truncated series of the wave function) set to 500 and L (the half- width of the well) set to 2. 
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Figure 5: The wave function for ground (left) and first excited (right) state for the bounded Morse 
potential, with depth parameter Vq = 400 and range parameter A = 1, where the results for the bounded 
case are obtained for / (the number of nonvanishing terms in the truncated series of the wave function) 
set to 500 and L (the half-width of the well) set to 2 
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4 Conclusion 

In this paper, we have applied the method of power-series expansion to a variety of 
symmetric one-dimensional anharmonic potential functions bounded by two infinite walls. 
We have compared our calculations of the low-lying energy levels with results from the 
inner-product quantization procedure. For all the potentials considered, we have chosen 
values for the number n of terms in the power series and the distance 2L between the 
walls that yield results of an accuracy superior to the results from the inner-product 
quantization procedure. More importantly, in analyzing the ground and first excited 
states of the deep double well, we have found the splitting of the energy levels to become 
apparent only after 26 significant digits for the mass parameter /i 2 = 25, in agreement with 
the result obtained from the inner-product quantization procedure. We have extended our 
calculation to the case of fx 2 = 35 and fx 2 = 50, where the splitting shows up after 42 and 
72 significant digits, respectively. For sextic, octic, dectic and duodectic oscillators, we 
have calculated the energy and wave function for the first excited states as well as for the 



ground states. Finally, we have applied our method to the well-known Morse potential, 
finding that our numerical results match the exact results. 

A reason for the capability of the present method is that for a bound state, the 
wave function is spatially localized. This means that the probability density |\l/| 2 has 
appreciable values in a finite region of space, outside which the probability density tends 
rapidly to zero. To obtain a good approximation, it is thus reasonable, as is shown by 
our numerical results for our chosen example of the Morse potential, to consider the 
corresponding problem in a finite interval, with a suitable width, bounded by two infinite 
walls. 

Although the two methods are equally powerful for the potentials considered in the 
present work, the power-series expansion in the finite range is convergent for potentials 
admitting a power-series expansion in the same finite interval. The power-series method 
has been also applied to the case of one dimensional multi-well oscillator [12]. In addition, 
it has been used in [13] to justify numerical results based on the Zinn- Justin conjecture 
[14]. We show in work under preparation that the power-series method can be extended 
to the case of three-dimensional spherically symmetric potentials. 
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